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Using Dynamical Mean Field Theory (DMFT) and Monte Carlo (MC) simulations, we study the 
ferromagnetic transition temperature (T c ) of a two-band model for Diluted Magnetic Semiconductors 
(DMS), varying coupling constants, hopping parameters, and carrier densities. We found that T c is 
optimized at all fillings p when both impurity bands (IB) fully overlap in the same energy range, 
namely when the exchange couplings J and bandwidths are identical. The optimal T c is found 
to be about twice larger than the maximum value obtained in the one-band model, showing the 
importance of multiband descriptions of DMS at intermediate J's. 

PACS numbers: 71.27.+a, 75.50.Pp, 75.40.Mg 



I. INTRODUCTION 

DMS are semiconductors where a fraction x(~0.0l- 
0.1) of nonmagnetic elements is replaced by magnetic Mil 
ions. The class of III-V DMS has recently attracted con- 
siderable attention after the experimental observation of 
high T c 's, due to significant improvements in molecular 
beam epitaxy techniques^ These compounds can play an 
important role in spintronic devices, and they represent 
a challenge to theory due to the combined presence of 
correlations and disorder. 

Using a one-band model for DMS, (i) the weak- 
coupling quadratic dependence of T c with J is captured 
correctly by DMFT—, which also provides bulk limit re- 
sults, while (ii) the MC techniques on finite clusters prop- 
erly handle the random Mn distribution and unveiled the 
reduction of T c at large J due to carrier localization— For 
J comparable to the hopping t, both techniques reached 
similar conclusions, showing that they can complement 
towards a comprehensive understanding of DMS. How- 
ever, much theoretical work remains to be done, particu- 
larly the consideration of the several bands active in real 
DMS materials beyond mean field approximations— 

In this paper, we study a two-band model for DMS 
at realistic low dopings x<^l, using a combination of 
nonperturbative techniques that gives unique character- 
istics to our work. To properly analyze the intermedi- 
ate and strong coupling J regime, MC simulations that 
correctly handle the random Mn distribution are cru- 
cial. Here, MC results for DMS models including more 
than one band are presented for the first time— On the 
DMFT side, pioneering calculations for two-band double- 
exchange Hamiltonians already exist— However, those 
calculations focused on special cases, while our current 
DMFT effort is more general, with two s=l/2 bands and 
arbitrary couplings, hoppings, and carrier densities. One 
of our main results is that T c can be substantially raised 
by considering multiband systems since at intermediate 
couplings the maximal T c at carrier filling p=x is ap- 
proximately twice larger than the highest T c obtained in 
the single-band model at filling p=x/2 (note that the al- 



ternative notation ph—p/x, where ph is the hole density 
as a fraction of p, will be used in some portions of the 
paper). The qualitative reason is that the IB cooperate 
to raise T c for values of the chemical potential (i where 
those bands partially or fully overlap— 

The paper is organized as follows: Section II describes 
the model; the DMFT results are presented in Section 
III, while Section IV is devoted to the Monte Carlo study. 
Conclusions and final remarks arc in Section V. 

II. MODEL 

It is known that in Mn-doped GaAs, the Mn ions sub- 
stitute for Ga cations and contribute itinerant holes to 
the valence band. The Mn ions have a half-filled d-shell 
which acts as a S = 5/2 local moment. Due to a strong 
spin-orbit (SO) interaction, the angular momentum L of 
the p-like valence bands mixes with the hole spin degree 
of freedom s and produces low- and high-energy bands 
with angular momentum j = 1/2 and 3/2, respectively. 
A robust SO split between these bands causes the holes 
to populate the j = 3/2 state, which itself is split by the 
crystal field into a rrij = ±3/2 band with heavy holes 
and a rrij = ±1/2 band with light holes. This is the rea- 
son why we choose to study two bands since this is the 
relevant number of orbitals in most III-V DMS. Since we 
do not work in a (j, rrij ) basis our Hamiltonian does not 
capture the orbital mixing in the Hund term— However, 
we roughly consider the diagonal SO effects in the mag- 
netic interactions by allowing different values of J in the 
two orbitals considered. The simple two-band model for 
DMS used here is given by the Hamiltonian 

W =-]L ^(Citi,a C ',J> + H - C -)-X] JlSl ' B V> W 

l,ij,a 1,1 

where 1=1,2 is the band index (not to be confused with 
angular momentum), i,j label sites (nearest neighbors 
for the hopping term), c^i a creates a hole at site i in 
the band I, s M = J2 a ,p c ti,a a aP c hi,P is tne s P m " op- 
erator of the mobile hole (a = Pauli vector), a and (3 
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FIG. 1: (Color online) T c versus the carrier concentration p, 
at various J2/J1, obtained with the DMFT technique. Here, 
x = 0.05, W1/W1 = 1, and J1/W1 = 0.5. The inset shows 
the corresponding DOS at T = 0. 



are spin indices, J; is the coupling between the core spin 
and the electrons of band I, and S/ is the spin of the 
localized Mn ion at randomly selected sites /, assumed 
classical in the MC simulations, ti is the hopping term in 
band I. The inter-band hopping ti2 (= £21) is zero at the 
nearest neighbor level in cubic lattices^. Even if ti2 is 
explicitly added, the conclusions arc similar as reported 
herein. While real DMS materials have zincblendc (ZB) 
structures, in this first nonpcrturbative study of a multi- 



band DMS model the simplicity of a cubic lattice allows 
us to focus on the dominant qualitative tendencies, a first 
step toward future quantitative studies with realistic ZB 
lattices. The model will be studied using DMFT and MC 
techniques. 



III. DMFT RESULTS 



A. Formalism and tests 

Within DMFT the self-energy is assumed to be local, 
S(p, iuj n ) — > S(?w„), assumption valid exactly only in in- 
finite dimensions. The information about the hopping of 
carriers on and off lattice sites, which are magnetic (with 
probability x) or nonmagnetic (with probability 1-x), is 
in the bare Green's function (GF) Qo(iu> n )^. The full 
GF G(iu) n ) can be solved by integration with the result: 
(G(iu n )) = x{[Go l {iu n ) + JSmo-]- 1 ) + (1 - x)(g (iu n )), 
where w n =(2n + \)ttT arc Matsubara frequencies. The 
average (X(m))= J dft m X(m)P(m) is over the local mo- 
ment m orientations, with probability P(m). The equa- 
tion above was derived for each band, and the equa- 
tions set for the GF were solved on a Bethe lattice 
with a bare semicircular nonintcracting density of states 
DOS=8(PF, 2 /4-e 2 ) 1 / 2 /7rVF ; 2 . In this case, (g^, 1 (««;„)) = 
z n - (W; 2 /16) (Gi(iu n )), where z n =iu n + fi, and Wi=4ti. 
T c is obtained by linearizing the bare inverse GF with 
respect to the local order parameter M = {m z ) m . To 
first order in M, the implicit equation for T c is: 
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where Bi(iui n ) satisfies a fourth-order-degree equation: completely separated from the valence bands. 



R3_ 

'- D Z 



(J 2 -I¥ 2 /I6)i? 2 +z n J 2 i? ( -(I-x)J 2 ^ 2 /4 



Oilil Our formula was tested in different cases: (1) for 
J 2 =0, we reproduce the one-band results of Ref0, and 
(2) at x=l, J2=0, and Ji — > oo we reproduced the results 
of Ref. [ill T c , contained in the Matsubara frequencies, 
is obtained from Eq.@ numerically. The equations for 
Bi, rewritten in real frequency space via iuj n — > ui, give 
the interacting DOS for each band at zero temperature 
when fi = 0H The solutions of these equations depend 
crucially on the ratio J;/ Wi . The critical value for the for- 
mation of well-defined IB, corresponding to carrier spins 
locally parallel to Mn spins, is Ji/Wi ~ 0.33. Below, the 
domain Ji/Wi < 0.33 is referred to as 'weak coupling', 
while 0.33 < Ji/Wi < 0.5 is the 'intermediate coupling'. 
The most interesting physics is observed at the bound- 
ary between these two regimes, i.e. when the IB are not 



B. DMFT critical temperatures varying exchange 
couplings 

In Fig.^ wc show T c vs. p, for different ratios J2/J1 
and at fixed W2/W\ = l and J\/W\ = 0.5, situation cor- 
responding to the existence of a well-defined 1=1 IB (al- 
though p<x in real DMS, in this paper the case P>x will 
also be studied for completeness, as done in RefQ). The 
inset shows the total interacting DOS evolution. The 
IB overlap if \ J 2 /W 2 - Ji/Wi| < 0.5. If the IB do not 
overlap, then each one determines T c separately, causing 
the double-peak structure observed for some J2/J1 ra- 
tios. The band with the largest Ji/Wi is filled first, for 
smaller //'s. At all p's, we found that T c is maximum 



3 



0.010 




0.025 0.050 0.075 0.100 0.125 0.15 
P 



FIG. 2: (Color online) Results obtained with the DMFT 
approximation: (a) T c vs. J2/J1, at W1/W2 — 1, and 
J1/W1 = 0.5, for the values of p indicated in (b). At p = 0.02 
and J2/J1 = 0, a finite T c /Wi ~ 0.0037 is caused by the 1=1 
band. At p = 0.05, T c is not zero for J2/J1 G (0.6,1.35), 
and it increases significantly due to band overlap. The case 
p — 0.07 (> x) corresponds to the first IB completely filled, 
(b) T c vs. J1/W1, at W1/W2 = 1 and J1/W1 = J 2 /W 2 , for 
the p's indicated, (c) T c vs. p at different ratios W2/W1, 
fixing Ji = J2 = 2. In all the frames x = 0.05. 



when J2/J1 = 1, namely when the IB fully overlap. The 
dependence of T c on J2/J1 at fixed p is in Fig.^a). The 
peak value is achieved when the bands fully overlap (i.e. 
at J2/Ji=l). Once the bands decouple, the value for T c 
matches one-band model results. 



C. DMFT critical temperatures varying 
bandwidths 

Let us consider now how changes in bandwidths influ- 
ence T c . In Fig. OJc), we show T c vs. p parametric with 
W2/W1, at fixed J1/W1 = 0.5 (intermediate coupling), 
and with J2/J1 = lM At small W2/W1 the second IB 
shall be located in a region of 10 smaller (i.e. farther 
from the valence bands), than the energy interval occu- 
pied by the I = 1 IB. Hence, the 1=2 IB will be the first to 
be filled. Decreasing ^2/^2, the second band moves to 
the right on the u> axis, towards the location of the first 
band. While the bands are still separated, each gives its 
own contribution to T c . The curves with W2/W1 = 0.35 
and 0.5 correspond to decoupled IB, while those with 
W2/W1 = 0.625 and 0.75 correspond to partially overlap- 
ping bands. Again, T c is maximized at all fillings when 
the bands fully overlap (W2/Wi = l), in good agreement 
with RefH Although W 2 /W 1 =l is not realistic in DMS, 
Mn doping materials with a relatively small heavy-light 
mass ratio will favor a higher T c . 



FIG. 3: (Color online) T c versus p at different ratios J1/W1, 
obtained using DMFT. The parameters x, Wi, W2, and Ji 
are as in Fig. 1. The inset contains the DOS at T = 0. 



D. Critical temperatures varying the exchange to 
bandwidth ratio 

Once established that T c is maximal for all p when 
^2/^1=1 an d W2/W\=l 1 let us analyze T c vs. p when 
J\jW\ varies. The results for T c , and total interacting 
DOS, are in FigCU At small coupling J x /W\ < 0.33, T c 
is small, flat, and much extended on the p axis, qual- 
itatively similar to the one-band results^ However, at 
intermediate coupling, T c is nonzero in the range from 
p = to p = 2x, adopting a parabolic form with the 
maximum at p = x, in contrast with the one-band model 
which gives a null T c when p = x. The explanation is 
straightforward: at p = x in the one-band model the IB 
is fully occupied leading to a vanishing T c , but for the 
same p in the two-band model both bands are half filled, 
which ultimately leads to the highest value for T c . The 
T c dependence on the ratio J\jW\ at some fillings p is 
displayed in Figl^b). At small coupling, T c correctly 
increases quadratically with J\/W\, but at strong cou- 
pling T c incorrectly continues growing, result which will 
be improved upon by the MC simulations shown below. 



IV. MONTE CARLO RESULTS 

The Hamiltonian Eq. Q was also studied numerically 
using MC techniques similar to those applied to Mn-oxide 
investigationsiiii^ The fermionic sector is treated exactly, 
while a MC simulation is applied to the classical localized 
spins. Cubic lattices with 5 3 and 6 3 sites were investi- 
gated. These lattice sizes have been shown to be suffi- 
cient for the comparison with DMFT results and to unveil 
general trends. In addition, the figures show only small 
variations for the T c estimations using the two lattices. 
However, if more sophisticated quantitative analysis is 
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FIG. 4: (Color online) (a) MC magnetization (in absolute 
value) vs. temperature (T), with J± — J2=4 and ti=t2=l, 
for the hole densities indicated. The dashed line is the ex- 
act asymptotic high-temperature value , which tends to 
only in the bulk limit. In this paper, the T c on the 5 3 cluster 
was (arbitrarily) defined as the T where M reaches ~5% of 
the 1 — Moo value (indicated by arrows in (a) and (b)). Other 
criteria lead to very similar conclusions, (b) Magnetization vs. 
T for Ph~i, ti=t2, and several (equal) magnetic couplings J; 
(c) Curie temperature vs. hole density for J=4 and different 
ratios of the band hoppings; (d) Curie temperature versus the 
(equal) magnetic interactions J for several hole densities and 
equal band hoppings ti and £2. Results for 5 3 (6 3 ) lattices 
are indicated by open (filled) symbols. Error bars due to the 
disorder (up to 7 samples) are only shown for a few points for 
clarity. 



required, clearly larger systems will be needed. One may 
suspect that actually the number of Mn spins may regu- 
late the size effects, rather than the number of sites. For 
the small values of x used in our study, the number of 
Mn spins is also very small and serious size effects could 
be expected. However, in practice this does not seem to 
occur, and moreover the comparison with DMFT shows 
similar results using both techniques. Perhaps in the 
small J's regime, the delocalized nature of the carriers 
smears the effects caused by the actual location of the 
Mn spins. These issues deserve further study, but for our 
purposes of unveiling trends in the multi-parameter space 
of DMS materials, the lattices here used are sufficient. 

Returning to the numerical data, the spin magnetiza- 
tion is the order parameter that was used to detect the 
ferromagnetic transition^ All quantities are in units of 
ti=l, and the density of magnetic impurities is x w 0.065. 
In Fig. 4(a), typical magnetization curves are presented 
at several carrier densities ph, and for Ji=J2=J=4 and 
£2=1. In excellent agreement with DMFT, it was ob- 
served that the estimated T c is the highest for p/,=l. The 
value of J used maximizes the critical temperature, and 
it was confirmed that it corresponds to the case where 
the IB are about to become separated from the valence 
band. In Fig. 4(b), it is shown how T c increases with J 
up to J=4, in agreement with DMFT. The strong cou- 
pling behavior is nevertheless different since T c decreases 



FIG. 5: (Color online) (a) T c vs. pn obtained with MC on a 
5 3 (6 3 ) lattice with £i=£2 = l for the values of J2/J1 indicated 
by the open (filled) symbols. Ji is fixed to 4, i.e. when the 
IB are about to separate from the valence band for band 1 
(inset Fig. 5(b)). For J 2 < Ji (e.g. J 2 /Ji=0.4 curve), T c is 
regulated by the IB in band 1, since the IB for band 2 is deep 
into the valence band. In this case, a single-band behavior is 
observed, with T c maximized for ph~0.5. For Ji > J\, both 
IB play a role. For J2 S> Ji (see J2/Ji=2), the two IB do not 
overlap, and for < Ph < 1 T c is determined by the band-2 
IB reaching a maximum at ph=0.5 and vanishing at ph~i- 
For larger ph, T c is controlled now by the IB 1, and it raises 
again passing through a maximum at ph=1.5 and vanishing 
at ph=2. For J2/Ji = 1.25, the two IB overlap and we observe 
residual local maxima at Ph=0.5 and 1.5, related to the single 
band physics, and a new local maximum at ph=l due to IB 
overlap for the corresponding value of the chemical potential. 
T c at ph=0.5 is boosted by the partial IB overlap as well, (b) 
T c vs. ph for £2/£i = l and several Js. Inset: low-T DOS. 



at large J. This is caused by hole localization in strong 
coupling^ beyond the capabilities of DMFT. The depen- 
dence of T c on the ratio of band hoppings is in Fig. 4(c), 
varying ph- These results are again in good qualitative 
agreement with DMFT (Fig. 2(c)). The maximum T c for 
all values of ph occurs when t<ilt\=\. However, when £2 
is very different from t\ the development of magnetism is 
regulated by only one of the IB and the results arc simi- 
lar to those obtained with a single band model, as shown 
in the curves for £2=0.35 and 0.625 in Fig. 4(c). For 
£2/£i closer to 1, a partial overlap of the IB occurs and a 
hump in T c develops at ph=l (see curve for £2=0.75). In 
Fig. 4(d) we show that, once £2/£i is optimized, a sim- 
ilar finite J maximizes T c for several phS. In all cases, 
the optimal J best balances the weak coupling behavior, 
with mobile holes not much affected by the interaction 
with the Mn ions, and the strong coupling region where 
the hole spins strongly align with the Mn spins, becom- 
ing localized. This "sweet spot" is achieved when the IB 
are about to be separated from the valence bands. 

T c vs. ph, at several ratios J2/J1 and for £i=£2 is pre- 
sented in Fig. 5(a). In agreement with DMFT (Fig. 1), 
T c is maximized for all values of ph if ^2=^1, with the 
highest value at f>/,=l. Overall, there is an excellent 
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agreement with DMFT, as described in the caption J^i 
T c vs. ph for several J\=J<i=J is in Fig. 5(b). At small 
J, again the MC results resemble those obtained with 
DMFT (Fig. 3). For, e.g., J=2 the IB are not formed 
yet (inset). In this regime, T c remains finite, although 
small, even for ph larger than 2. Increasing J, a nonzero 
T c is obtained only for ph between and 2, due to IB 
formation. T c reaches a maximum at J=4. 

V. CONCLUSIONS 

We have carried out the first study of a multiband 
model for DMS using a powerful combination of non- 
perturbative techniques, DMFT and MC. We found the 
parameter regime that maximizes T c . This happens at 
intermediate couplings and for all hole densities when 
J\l Ji=\ and VFi/W^l. The maximum T c is obtained 
at p = x, in contrast with the one-band model which has 
a vanishing T c at the same doping. In addition, T c at 
filling p = x/2 in the one-band case is smaller than with 
two bands by a factor ^2. In view of the simplicity of the 
main results, it is clear that adding an extra band to the 
calculations (which is relevant for system with negligible 
SO, but considerably raises the CPU cost) will only lead 
to a further increase in T c when all the IB overlap. 

The excellent agreement DMFT-MC is somewhat sur- 
prising due to the fact that Monte Carlo considers the 
influence of the random location of the Mn sites much 
better than DMFT. However, at small and intermedi- 
ate J's, the carriers can be sufficiently delocalized that a 
smearing effect may occurs and considering the quenched 
disorder only in average appears to be sufficient. Cer- 



tainly at large J's the MC and DMFT methods give to- 
tally different answers, with MC capturing the correct 
localization result. 

The approach described here is also quantitative. In 
fact, using GaAs realistic parameters such as p=0.005 
(p/^O.l/Mn), a bandwidth - 10 cV 2.5 eV), 

t2=(l/9)<i, and assuming J\jt\=\ and J\=J^ we ob- 
tain T c s=s 175 K, i.e. within the experimental range. 
While this excellent agreement with experiment may 
be accidental, the trends are reliable and the result im- 
proves upon single-band estimations. Moreover, for opti- 
mal t2/ti=l and Ji/<i=2, the T c raises to ~ 340 K, even 
at small ph=0.1, setting the upper bound for DMS under 
a two-band model description using a cubic lattice. 

The general qualitative picture presented here can be 
used to search for DMS with even higher T c 's than cur- 
rently known. Our results suggest that semiconductors 
with the smallest heavy to light hole mass ratio, such as 
AlAs, could have the highest T c if the couplings J could 
be tuned to its optimal value. The present effort paves 
the way toward future nonperturbative studies of DMS 
models using realistic ZB lattices, and points toward 
procedures to further increase the Curie temperatures. 

VI. ACKNOWLEDGMENTS 

We acknowledge conversations with K. Aryanpour, M. 
Berciu, R.S. Fishman, J. Moreno, and J. Schlicmann. 
This research is supported by grant NSF-DMR-0443144. 
ORNL is managed by UT-Battelle, LLC under Contract 
No. De-AC05-00OR22725. MC simulations used the 
SPF program ( http :/ /mri- fre.ornl.gov /spq ) . 



1 H. Ohno, Science 281, 951 (1998); T. Dietl, H. Ohno, F. 
Matsukura, J. Cibert, and D. Ferrand, Science, 287, 1019 

(2000) ; S. J. Potashnik, K.C. Ku, S.H. Chun, J.J. Berry, 
N. Samarth, and P. Schiffer, Appl. Phys. Lett. 79, 1495 

(2001) . 

2 A. Chattopadhyay, S. Das Sarma, and A.J. Millis, 
Phys. Rev. Lett. 87, 227202 (2001), and refs. therein. For 
a DMFT review, see A. Georges, G. Kotliar, W. Krauth, 
and M.J. Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 

3 G Alvarez, M. Mayr, and E. Dagotto, Phys. Rev. Lett. 
89, 277202 (2002); ibid, Phys. Rev. B 68, 045202 (2003). 

4 J. Sinova, T. Jungwirth, S.-R. Eric Yang, J. Kucera, and A. 
H. MacDonald, Phys. Rev. B66, 041202 (2002); T. Jung- 
wirth, J. Konig, J. Sinova, J. Kucera, and A.H. MacDonald 
, ibid 66, 012402 (2002). 

5 In J. Schliemann, J. Konig and A.H. MacDonald, Phys. 
Rev. B64, 165201 (2001), MC studies of multiband models 
at small J were presented, but without using an underlying 
lattice, and replacing the fermion trace by a ground state 
expectation value. 

6 K. Aryanpour, J. Moreno, M. Jarrell, and R.S. Fishman, 
Phys. Rev. B7 2 045343 (2005); J. Moreno, R.S. Fishman, 
and M. Jarrell, |cond-mat/0507487| 



7 The term "impurity band" is widely used in the DMS liter- 
ature (see e.g. M. Berciu and R. N. Bhatt, Phys. Rev. Lett. 
87, 107203 (2001); J. Okabayashi, A. Kimura, O. Rader, 
T. Mizokawa, A. Fujimori, T. Hayashi, and M. Tanaka, 
Phys. Rev. B64, 125304 (2001)), although it should not be 
confused with the same term used in nonmagnetic doped 
semiconductor literature. Perhaps "exchange split valence 
band" could be a better terminology. 

8 G Zarand and B. Janko, Phys. Rev. Lett. 89, 047201 

(2002) . In Ref0 it was shown that the frustration can be 
turned off by giving equal masses to the light and heavy 
sectors. Our emphasis has been on this particular case, 
since under this circumstance the T c is maximized. 

9 A. Moreo, to be published. 

Details will be published elsewhere. 

1 R.S. Fishman and M. Jarrell, Phys. Rev. B 67, 100403 

(2003) . 

2 The interacting DOS(o;) is -(2/tt) Yjf =1 lm[Bi(uj)]. 

3 The relative position of the bands is fixed by shifting nu- 
merically the 1—2 valence band such that both valence 
bands start at the same energy, to mimic the degeneracy 
of the light and heavy bands of GaAs at T. 

4 E. Dagotto, T. Hotta, and A. Moreo, Phys. Rep. 344, 1 



(2001). rescale the hopping in DMFT, is considered. 

This agreement is even quantitative for the optimal T c and 
J coupling, once the factor \f2d (d=dimension), use to 



